Magnetic heat transport in i? 2 Cu04 (R = La, Pr, Nd, Sm, Eu, and Gd) 
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We have studied the thermal conductivity k on single crystalline samples of the antiferromagnetic 
monolayer cuprates .R2CUO4 with R — La, Pr, Nd, Sm, Eu, and Gd. For a heat current within 
the CUO2 planes, i. e. for K a b we find high-temperature anomalies around 250 K in all samples. 
In contrast, the thermal conductivity k c perpendicular to the CuCh planes, which we measured 
for R = La, Pr, and Gd, shows a conventional temperature dependence as expected for a purely 
phononic thermal conductivity. This qualitative anisotropy of Ki and the anomalous temperature 
dependence of K a b give evidence for a significant magnetic contribution ft mag to the heat transport 
within the Cu02 planes. Our results suggest, that a large magnetic contribution to the heat current 
is a common feature of single-layer cuprates. We find that K mag is hardly affected by structural 
instabilities, whereas already weak charge carrier doping causes a strong suppression of ft mag . 



PACS numbers: PACS numbers: 74.72.-h, 66.70.+f 

I. INTRODUCTION 

Various low-dimensional spin systems show an unusual 
thermal conductivity k with a double-peak structure as 
a function of temperature. There is growing evidence 
that this anomalous behavior arises from magnetic ex- 
citations contributing to the heat transport. The most 
clear experimental evidence is found in the spin-ladder 
compounds, where a double-peak structure with a huge 
high-temperature maximum of k is present for a heat 
current parallel to the ladder direction, but absent for 
a heat current perpendicular to the laddersi^ For one- 
dimensional spin-chain systems the experimental results 
are less clear. For the spin-chain compounds SrCu02 
and Sr2CuC>3 Sologubenko et al. find a sizeable ex- 
tra contribution to k along the chain direction, which 
is missing in the other directions^ In the spin-Peierls 
compound CuGeOa k along the chain direction has two 
low-temperature maxima and one of them was attributed 
to a magnetic contribution in Ref. 4. However, this inter- 
pretation is questionable, because a similar double peak 
is also present in k perpendicular to the chain direction. 5 
For the Haldane-chain (S = 1) System AgVP2S6, a mag- 
netic contribution seems to play a role toOf& but the abso- 
lute values are much smaller than in the S = 1/2 systems. 
The results in ID systems raise the question whether 
a sizeable heat current due to magnetic excitations is 
also present in two-dimensional magnets. This was dis- 
cussed for the low-temperature thermal conductivity of 
K 2 V 3 8 and Nd 2 Cu0 4 i£&£ The latter is one of the in- 
sulating parent compounds of high-temperature super- 
conductors containing Cu02 planes, which represent the 
perhaps most studied two-dimensional antiferromagnets 



so farifi Whereas the studies on Nd 2 Cu0 4 (Refs. HI 
mainly concern the magnetism of the Nd 3+ moments, 
the influence of the Cu 2+ moments is present at higher 
temperature. In the layered perovskite La2Cu04, the 
thermal conductivity n a b for a heat current along the 
CuC>2 planes exhibits a pronounced double-peak struc- 
ture with a low-temperature maximum around 25 K and 
a second one around 250 K. In contrast, the thermal con- 
ductivity k c perpendicular to the CuC>2 planes has only 
one low-temperature peakiii*i2iiii4 These findings have 
been interpreted in terms of an additional heat trans- 
port parallel to the CuC>2 planes due to magnetic exci- 
tations. However, a double-peak structure can also be 
explained by phonons only. Any additional scattering 
mechanism which acts in a narrow temperature range 
suppresses k in that temperature window and as a result 
K may exhibit two peaks. For example, in SrCu2(B03)2 
a double-peak structure is caused by resonant scatter- 
ing of acoustic phonons by magnetic excitations^ Such 
a mechanism does not apply for La2CuC>4. However, 
La2Cu04 has a structural instability with low-lying opti- 
cal phonon branches, which could also serve as scatterers 
for the acoustic phonons. Such an explanation has been 
proposed by Cohn et al. for the heat transport data of 
YBaCu06+5, which show a similar temperature depen- 
dence of K a b as La2Cu04ii& 

Recently, we have measured K a b of Sr2Cu02Cl2 in or- 
der to investigate the possibility that resonant phonon 
scattering due to the structural instability may cause the 
double-peak structure in n a b of La2CuC>4. S^CuO^Cls 
and La2CuC>4 are almost isostructural, but Sr2Cu02Cl2 
has no structural instability. Nevertheless we also found 
a second high-temperature maximum in K a b and took this 
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as evidence for the second peak of the in-plane heat con- 
ductivity being caused by magnetic excitations^ Bas- 
ing on this observation we expect that a pronounced 
magnetic contribution to the thermal conductivity is a 
common feature of the layered cuprates. However, the 
high-temperature peak of n a b of La2CuC>4 is significantly 
larger and its low-temperature peak is much smaller than 
the corresponding peaks observed in S^CuC^C^. These 
quantitative differences could arise from the absence or 
presence of a structural instability and/or weak charge 
carrier doping in the different samples. In order to inves- 
tigate whether the high-temperature peak of K a b is in- 
deed an intrinsic feature of the antiferromagnetic Cu02 
planes and to get more insight about the influence of 
structural changes, we have studied the thermal conduc- 
tivity on single crystals of -R2CUO4 with different rare 
earths R, which realize different structures. Some of 
these compounds even show a structural transition as a 
function of temperature pl&iiiSSiSi Up to now only the in- 
plane thermal conductivity for R = Pr and Nd has been 
studied; 8 i 9 i 22 i 23 and the studies of K a b and k c of La 2 Cu0 4 
have concentrated on the influence of doping»iiii2ii2ii4 

Here, we present measurements of K a b for R = Pr, Nd, 
Sm, Eu, and Gd and of k c for R = Pr and Gd. Our study 
clearly shows that the anomalous double-peak structure 
of K a b is present for all -R2CUO4 and confirms that a 
large magnetic contribution to the heat transport is an 
intrinsic property of the Cu0 2 planes. In addition, we 
find that this magnetic contribution is hardly affected by 
a structural instability while the phononic contribution 
is strongly suppressed. 



II. MAGNETIC AND STRUCTURAL 
PROPERTIES OF R 2 Cu0 4 

A common feature of La2Cu04 and -R2CUO4, with 
R = Pr, Nd, Sm, Eu, and Gd is the layered structure 
with planes consisting of a Cu02 square lattice. These 
planes are a good realization of a two-dimensional an- 
tiferromagnetic S = I Heisenberg square lattice. Bi- 
magnon Raman scattering 24 yields exchange constants 
J between w 1200 - 1400 K for T = 300 K given in 
TableHJ 25 Finite inter-plane couplings J± lead to three- 
dimensional antiferromagnetic ordering with Neel tem- 
peratures of Tn ~ 250 . . . 320 K (see Table[IJ|. In general, 
the ordering temperature is determined by the ratio of 
the inter- and the intra-plane coupling. However, crystal 
quality and the oxygen stoichiometry strongly influence 
Tn, too. For example, very small amounts of excess oxy- 
gen drastically suppress Tn of L^CuCU+.s* 26 Apart from 
the magnetic Cu subsystem, the compounds with mag- 
netic rare earth ions R contain another magnetic sub- 
system. In all these compounds the behavior of the Cu 
subsystem is very similar. However, the details of the 
magnetic structure are determined by the competition 
between the different couplings (Cu-Cu, R-R, i?-Cu)i 27 i 28 

La 2 Cu0 4 crystallizes in the so-called T structure (also 



called K2MO4 structure). The CUO4 plaquettes of the 
planes and the apex oxygen ions form CuOg octahe- 
dra. At high temperatures La2Cu04 is in the high- 
temperature tetragonal phase (HTT phase). At 530 K 
a structural phase transition takes placed where the oc- 
tahedra tilt leads to the low-temperature orthorhombic 
phase (LTO), which is stable down to lowest temperature. 
Due to the octahedron tilt the point bisecting the nearest 
neighbor Cu-Cu distance is no longer a center of inversion 
symmetry giving rise to a Dzyaloshinski-Moriya (DM) 
type interaction. Si^CuO^C^ is almost isostructural to 
La2Cu0 4 with the La 3+ ions being substituted by Sr 2+ 
and the apex 2 ~ by CI" ions. In Sr 2 Cu0 2 Cl2 the HTT 
phase is stable down to the lowest temperature and due 
to inversion symmetry no DM exchange and consequently 
no spin canting occurs. The -R2CUO4 compounds crystal- 
lize in the tetragonal so-called T'-structure. While the T 
structure may be viewed as a stacking of one Cu02 layer 
followed by two LaO (or SrCl) layers, the stacking of 
the T' structure is one CuC>2 layer followed by a layer of 
R 3+ ions, a layer of 2 ~ ions, and finally another layer 
of R 3+ ions. Consequently, there are no apex oxygen 
ions present in the T' structure and the basic building 
blocks are CUO4 plaquettes instead of the CuOe octahe- 
dra of the T structure. For R = Pr, Nd, and Sm the 
T' structure is stable over the entire temperature range. 
For EU2CUO4 and Gd2Cu04 structural phase transitions 
are observed with transition temperatures of 170 K and 
685 K, respectivelyi22i2i The structural changes can be 
described by an alternating rotation of the CUO4 pla- 
quettes around the c axis. The rotation angles amount 
to 2.3° at 20 K for R = Eu and to 5.2° at 300 K for 
R = Gdi22i2i The structural transitions transform both 
the T and the T' structure into orthorhombic structures, 
but all these crystals are usually strongly twinned with 
respect to the a and the b axes. As in La2Cu04, the 
lower symmetry in the distorted T' phases gives rise to 
a DM interaction and leads to a canting of the magnetic 
moments. 



III. EXPERIMENTAL 

The -R2CUO4 single crystals were grown in Pt cru- 
cibles by the top-seeded solution method as described in 
Ref-IU The La2CuC>4 crystal was grown by a traveling- 
solvent floating zone method. The finite DM interaction 
for R — La and Gd causes a weak ferromagnetic mo- 
ment, which allows an easy determination of the Neel 
temperature Tn by measurements of the magnetic sus- 
ceptibility x- F° r R = Gd we find Tn — 292 K in agree- 
ment with the highest values reported for Gd2Cu04i^ 3 ^ 
Our La2Cu04+5 single crystal has a comparatively low 
Tn = 245 K due to excess oxygen. A comparison of 
the in-plane resistivity (not shown) with Ref. |3^ yields 
S w 0.01. The determination of Tn from x(T) does nei- 
ther work in the undistorted samples with R = Pr, Nd, 
and Sm, since there is no weak ferromagnetism nor for 
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TABLE I: Exchange constants J at 300 K (Ref. 1241) . Debye temperatures Od, sound velocities v s , sample sizes, and Neel temperatures Tn of 
H2CUO4 (see text and the respective references). P, D, U, and u arc fit parameters for the fits of the phononic contribution of k. P an d D 
describe the scattering on point defects and planar defects, respectively, whereas U and u model Umklapp scattering (for details see Ref. El). If 
available, the values for Qjj and v s are taken from literature, otherwise similar values have been used for the fits. 
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R = Eu, where the structural transition takes place well 
below Tn . According to Ref. there is a slope change 
of Xab at Tn in NCI2CUO4, but we could not reproduce 
such a feature in \ab of our crystal. For Pr 2 Cu04 we 
determined Tn — 250 K at the Laboratoire Leon Bril- 
louin, Saclay by neutron diffraction, which is well below 
the maximum values up to Tn — 280 K reported for this 
compoundiffiSi^ 9 . Unfortunately, our Nd2CuC>4 crystal 
is too small for neutron diffraction and this method can- 
not be applied for R = Eu and Sm because of the large 
neutron absorption cross section of these elements. 

All crystals have been oriented using a Laue camera 
and cut into rectangular pieces. Sample sizes are listed in 
Table [I] The accuracy of the orientation with respect to 
the crystal axes is about 2°. Since the Sm 2 Cu04 crystal 
has approximately the shape of a cuboid, it has not been 
cut. Here, the misalignment amounts to ~ 10°. The 
shape of Pr2Cu04 and Gd2CuC>4 allowed measurements 
of k c with a heat current ja parallel to the c axis, i. c. 
perpendicular to the CuC>2 planes, and of n a b with jn 
within the CuC>2 planes. For R = Nd, Sm, and Eu, we 
could only measure K a0l because these crystals were very 
thin with lengths of less 0.4 mm parallel to the c axis. For 
R = Pr, Nd, and Gd we measured K ao with jn parallel 
to the a! axis of the HTT phase, which has an angle of 
45° with respect to the orthorhombic a and b axes. For 
R = Eu and Sm ju had an arbitrary orientation with 
respect to the a and b axes. 

The thermal conductivity has been measured by a 
standard steady-state method. One end of the sam- 
ple has been attached to the sample holder by silver 
paint and a small resistor has been glued to the opposite 
end of the sample by an insulating varnish (VGE-7031, 
LakeShore). The temperature of the sample holder has 
been stabilized and an electrical current through the re- 
sistor has been used to produce a heat current through 



the sample. The resulting temperature gradient has 
been determined by a differential Chromel-Au+0.07%Fe- 
thermocouple, which has been also glued to the sam- 
ple. Typical temperature gradients were about 0.5% of 
the sample temperature. The absolute accuracy of our 
method is restricted to about 10% because of uncertain- 
ties in determining the sample geometry and, in particu- 
lar, the exact distance between the two ends of the ther- 
mocouple. The relative accuracy is about one order of 
magnitude better. Radiation losses are negligible for our 
sample geometry. 



IV. RESULTS 

A. Gd 2 Cu0 4 and Pr 2 Cu0 4 

FigureQ] shows the thermal conductivity of Gd2Cu04 
measured on two different samples. One sample allowed 
to measure K a b and re c , whereas the second crystal was 
too thin to measure n c . In Figure[5]we display K a i, and 
n c of Pi'2Cu04. Apart from differences around the low- 
temperature maximum, which is very sensitive to the 
sample quality, our data agree well to the previously re- 
ported K a b of Pr 2 Cu04i2 3 . For both compounds k c follows 
the typical temperature dependence of the thermal con- 
ductivity of acoustic phonons. As shown by the solid 
lines, k c can be reasonably well described within a De- 
bye model. We use the same model and nomenclature 
as in our previous publication (see Eqs. (1) and (2) of 
Ref. Il7|) and the corresponding fit parameters are given 
in Table [I] In both crystals K a b exceeds k c over the 
entire temperature range. The low-temperature max- 
ima for Gd2Cu04 are slightly shifted in temperature. 
The most striking anisotropy is, however, the additional, 
broad maximum of n ao around 250 K. Although the low- 
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FIG. 1: (Color online) In-plane (« c ) and out-of-plane (K a b) 
thermal conductivity of Gd2Cu04. K a b was measured on two 
different crystals. Solid lines are fits by the Debye model and 
the dashed line is k c multiplied by a factor of 1.25 (see text). 



FIG. 3: (Color online) In-plane thermal conductivity of 
Nd 2 Cu0 4 , Sm 2 Cu0 4 , and Eu 2 Cu0 4 . For R = Eu the low- 
temperature maximum is suppressed. The inset compares K a b 
(•) of our Nd 2 Cu0 4 crystal with data (o) from Ref.[g. 



temperature maxima of K a b strongly differ for the two 
crystals with R = Gd indicating differences in the crystal 
quality, the magnitudes of their high-temperature max- 
ima are almost identical. In Pr2Cu04 n a b also shows 
an additional high-temperature maximum, but its mag- 
nitude is less pronounced (see section Ivl) . 

The double-peak structures of n a b cannot be modeled 
by the usual Debye model, but it is possible to describe 
the low-temperature maxima up to about 50 K. Com- 
paring the corresponding fit parameters of n a b and k c 
(see Table [J, the largest differences are found for the pa- 
rameter D, which is significantly larger for the fits of k c 



than for those of n a b- This parameter gives the strength 
of phonon scattering by planar defects and it appears 
reasonable that due to the layered structure of .R2CUO4 
scattering by planar defects should be more effective for 
a heat current perpendicular to the planes than for ju 
within the planes. Thus we interpret the different mag- 
nitudes of the low-temperature maxima of k c and n a b 
as a consequence of the layered structure. The high- 
temperature maxima of K a b around 250 K will be dis- 
cussed in sectionlVl 
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FIG. 2: (Color online) In-plane (k c ) and out-of-plane (K a b) 
thermal conductivity of Pr 2 Cu04. Lines are fits by the Debye 
model (see text). Inset: The same data on double-logarithmic 
scales with k c multiplied by a factor of 1.2. The temperature 
dependencies of n a b and k c around the low-temperature max- 
ima are nearly the same. Above about 70 K k c follows a 1/T 
behavior (solid line) , whereas an anomalous high-temperature 
upturn is present in K ab . 



B. Nd 2 Cu0 4 , Sm 2 Cu0 4 , and Eu 2 Cu0 4 

In Fig.|3| we show the in-plane thermal conductivities 
K a b of Nd2Cu0 4 , Sn^CuGj, and EU2CUO4. The data for 
R = Nd and Sm are very similar to each other and also 
to those of R = Pr and Gd. In all crystals n a b exhibits 
a well defined low-temperature peak around 20 K and an 
additional broad maximum around 250 K. For R = Nd 
similar results have been obtained by Jin et al&, but K a b 
of our crystal is systematically lower in the entire tem- 
perature range (see Inset of Fig.|3J). For Eu 2 Cu04 the 
low-temperature peak of n a b is almost completely sup- 
pressed. As described above, for R = Eu the structural 
transition takes place at 170 K. The suppression of the 
low-temperature peak is most probably a consequence 
of this structural instability, which prevents a strong in- 
crease of the phonon mean free path at low temperatures. 
In contrast, however, the high-temperature maximum of 
K a b is hardly affected by this structural transition. The 
high-temperature maximum for EU2CUO4 is even more 
pronounced than for the structurally stable crystals with 
R = Nd and Sm. 



5 



V. DISCUSSION 

Our data of k c of Pr 2 Cu 04 and Gd 2 Cu04 together 
with k c of La 2 Cu0 4 (Refs. lllll2ll3ll4T) clearly reveal 
that, on the one hand, the out-of-plane thermal conduc- 
tivities of all these compounds of slightly different struc- 
tures (tetragonal and orthorhombic T', and orthorhom- 
bic T) do not exhibit any indications of an anoma- 
lous high-tcmperature contribution. On the other hand, 
all in-plane conductivities clearly show additional broad 
high-temperature maxima. The fact that these high- 
temperature maxima are present in crystals without 
(i? 2 Cu0 4 with R = Pr, Nd, Sm, and Sr 2 Cu0 2 Cl 2 ) or 
with (different) structural instabilities (T: R = La; T': 
R = Eu and Gd) unambiguously shows that the anoma- 
lous contribution to K a b does not depend on the exis- 
tence of structural instabilities. This complements our 
previous suggestion based on a comparative study of k 
of Sr 2 Cu0 2 Cl 2 and La 2 Cu04 and, as we have discussed 
in detail in Ref. 0, the most natural explanation for 
the high-temperature maximum of K a b is an additional 
contribution to the heat transport caused by magnetic 
excitations. 

The similar behavior of k for all i? 2 Cu04 confirms that 
such a magnetic contribution n mag to the in-plane heat 
transport is indeed an intrinsic property of the Cu0 2 
planes. It remains, however, to clarify what determines 
the magnitude of K mag . For a quantitative analysis we 
consider the in-plane thermal conductivity as the sum of 
a phononic and a magnetic contribution 



^mag , 



(1) 



which are only weakly coupled to each other. In general, 
such an Ansatz can be used when the characteristic en- 
ergy scales for the two contributions are well separated 
from each other. For example, this is usually the case for 
electronic and phononic heat transport since the Fermi 
temperature is much larger than the Debye temperature, 
i.e. Tp 3> &d- I n the case of i? 2 Cu04 it is a priori not 
clear whether the assumption of weakly coupled phononic 
and magnetic contributions is fulfilled, since the magnetic 
coupling J is only about four times as large as 6c (see 
Table^l. However, the experimental observation that K a b 
shows two characteristic maxima, which are well sepa- 
rated from each other, encourages us to use Eq. QJ. 

In order to separate K mag from K a \, we assume that 
in the region of the low-temperature peak K mag is neg- 
ligibly small and fit the data for T < 50 K (Gd 2 Cu04: 
T < 85 K) by the Debye model (see Eq. (1) of Ref. E3) 
and subtract the extrapolation of the fit from the mea- 
sured data up to room temperature, i. e. K mag = n a b — 
«£k 4SI This analysis is not possible for Eu 2 Cu0 4 , because 
the low-temperature maximum is not well-enough pro- 
nounced. However, even without a fit we expect rather 
similar values of K mag for R = Nd, Sm, and Eu, be- 
cause the high-temperature data of K a b are very similar 
for these crystals (see Fig-El - One can check the ap- 
plicability of the Debye model to describe the phononic 
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FIG. 4: (Color online) Magnetic contributions to the in-plane 



thermal conductivity, calculated via k„ 



where 



K p h is determined by a Debye fit of the low-temperature maxi- 
mum. Upper panel: Values calculated from our measurements 
of Kab °f JfeCuCU and from the data from Ref. 0. Lower panel: 
The same analysis for various data of La2Cu04+^ taken from 
and of our crystal with Tn = 245 K. Inset: The 
maximum of the calculated K mag vs. the Neel temperature for 
La2Cu04+i (see text). 



contribution by corresponding fits of k c , i.e. we restrict 
the temperature range of the fit to the low-temperature 
maxima and then compare the high-temperature extrap- 
olations of the fits to the measured n c . As shown in 
Figs.n and |3 these fits yield a good description of n c 
for R = Gd over the entire temperature range, whereas 
in the case of Pr the high-temperature values of K/q are 
slightly underestimated by the fit. This probably arises 
from the sharper low-temperature peak for R = Pr. Since 
the low-temperature peaks of n a b for R = Pr, Nd, and 
Sm are also rather sharp, one may expect that the cor- 
responding Debye fits will also underestimate the high- 
temperature values of the phononic contribution of n a b 
and consequently the magnetic contributions K mag may 
be overestimated to some extent. 

In the upper panel of Fig.|H we compare the resulting 
K ma g of all our i?2Cu04 crystals, our previous result of 
Sr 2 Cu0 2 Cl 2 (Ref.Ujj), and re mag obtained from an anal- 
ysis of K a b measured on Nd 2 Cu04 (Ref.0). For compari- 
son, we also show K mag obtained from the literature data 
of various La 2 Cu04 crystalsiii*i2iiii Obviously, the tem- 
perature dependence of K ma g is very similar for all crys- 



6 



tals, but the magnitude of the broad maximum varies 
between about 7 and lOW/Km for our Sr2Cu0 2 Cl 2 and 
i?2CuC>4 crystals and from about 12 to 25W/Km for 
the various crystals from literature with R = La and 
NdAiiii2ii£ Although these differences are not too large, 
we do not think that they simply arise from the experi- 
mental uncertainty in the quantitative determination of 

^mag ■ 

Because the magnetic properties are rather similar for 
the different crystals, we expect that the different K mag 
mainly arise from differences in the scattering of the 
magnetic excitations. Possible scattering mechanisms 
are scattering between magnetic excitations and scat- 
tering by defects, phonons, and charge carriers. One 
may suspect that scattering between magnetic excita- 
tions, comparable e. g. to phonon-phonon Umklapp scat- 
tering, plays the most important role with respect to the 
temperature dependence of n mag . A deeper analysis of 
this scattering requires a detailed theoretical model for 
the dynamics of magnetic excitations, but even without 
such a model one may conclude that the similar magnetic 
properties naturally explain the similar temperature de- 
pendencies of K mag of the different compounds. 

The influence of defects and charge carriers on the 
thermal conductivity has been investigated in Zn- and 
Sr-doped La2Cu04i 13 i 14 . It has been found that Sr dop- 
ing suppresses fc mag much stronger than Zn doping and 
kmag vanishes almost completely above about 1% Sr. For 
both dopings the magnetic system is diluted, either by 
replacing magnetic Cu 2+ by nonmagnetic Zn 2+ ions or 
by the formation of Zhang- Rice singlets due to the intro- 
duced holes. However, the mobility of the holes strongly 
enhances the effect of charge-carrier doping, what is also 
reflected in a much stronger suppression of Tn by Sr dop- 
ing as compared to Zn doping4i As shown in the lower 
panel of Fig.01 the magnitude of K mag for the various 
nominally undoped La2CuC>4 crystals varies by about a 
factor of two. Since it is known that La2Cu04+a is likely 
to have some excess oxygen, we plot the maximum of 
Kmag as a function of Tn , which is very sensitive to small 
amounts of S (Inset of Fig.^J. The observed correla- 
tion between the magnitude of K mag and Tn is a clear 
indication that the different magnitudes of K mag arise 
from small amounts of charge carriers in the different 
crystals. 12 i 42 

One may suspect that the higher values of /-trna.p- of 
i? 2 Cu0 4 with R = La and Nd from Refs. I8lllll2ll3l 

in comparison to our crystals could result from a weak 
charge carrier doping in our crystals. The rather low 
Tn — 250 K of our Pr2Cu04 crystal compared to the 
T N values up to ~ 280 K reported in literaturo 37 i 38 i 39 
supports this view. However, this argumentation can 
neither explain the low K mag of our Gd2CuC>4 with a 
large T N ~ 292 K nor does it hold for Sr 2 Cu0 2 Cl2, 
which is commonly believed to be very stable with re- 
spect to charge carrier doping. Unfortunately, not much 
is known about possible variations of the oxygen content 
in i?2Cu04. Irrespective of the question of the exact oxy- 



gen stoichiometry, our finding that K mag is very similar in 
crystals with and without structural instabilities leads to 
the conclusion that scattering by phonons seems to play 
a minor role for the magnetic heat transport in the Cu02 
planes. This is most clearly seen in EU2CUO4 where the 
phononic low-temperature peak of EU2CUO4 is strongly 
suppressed by a structural instability whereas its mag- 
netic high-temperature maximum is hardly affected. 



VI. CONCLUSIONS 

In summary, we have studied the thermal conductivity 
of the rare earth cuprates -R2CUO4 for both, a heat cur- 
rent perpendicular (R = Pr and Gd) and parallel (i? = 
Pr, Nd, Sm, Eu, and Gd) to the Cu0 2 planes. The out- 
of-plane thermal conductivity shows the typical tempera- 
ture dependence of a purely phononic thermal conductiv- 
ity with a low-temperature maximum, whose magnitude 
depends on the crystal quality. In contrast, the in-plane 
conductivity for all crystal exhibits a pronounced double- 
peak structure consisting (i) of a low-temperature peak 
similar to that of the out-of-plane thermal conductivity 
and (ii) of an anomalous high-temperature contribution 
with a broad maximum around 250 K. Such an anisotropy 
between the in-plane and the out-of-plane thermal con- 
ductivity is also found in LazCuOV 1 ! 12 ! 13 ! 14 The fact 
that the double-peak is present in the structurally sta- 
ble -R2C11O4 with R = Pr, Nd, and Sm unambiguously 
rules out the possibility that the double-peak structure 
is caused by a structural instability which is present for 
R = La, Eu, and Gd. The qualitative anisotropy be- 
tween the in-plane and the out-of-plane thermal con- 
ductivity and the rather similar high-temperature be- 
havior of the out-of-plane thermal conductivity for all 
the different crystals gives clear evidence that this addi- 
tional high-temperature contribution arises from a size- 
able heat transport by magnetic excitations within the 
Cu02 planes. Our analysis yields a magnetic contribu- 
tion to the in-plane thermal conductivity between about 
7 to 25W/Km depending on the R system. In weakly 
doped La2Cu04 this magnetic contribution is strongly 
suppressed showing that scattering of magnetic excita- 
tions by mobile charge carriers plays an important role. 
In contrast, the structural instability does hardly influ- 
ence the magnetic thermal conductivity indicating that 
scattering of magnetic excitations by soft or anharmonic 
phonons plays a minor role. In order to clarify the role of 
scattering between magnetic excitations theoretical mod- 
els describing the dynamics of magnetic excitations would 
be highly desirable. 
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